Visualization of correlation cascade in spatio-temporal chaos using wavelets 
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We propose a simple method to visualize spatio-temporal 
correlation between scales using wavelets, and apply it to two 
typical spatio-temporally chaotic systems, namely to coupled 
complex Ginzburg-Landau oscillators with diffusive interac- 
tion, and those with non-local interaction. Reflecting the dif- 
ference between underlying dynamical processes, our method 
provides distinctive results for those two systems. Especially, 
for the non-locally interacting case where the system exhibits 
fractal amplitude patterns and power-law spectrum, it clearly 
visualizes the dynamical cascade process of spatio-temporal 
correlation between scales. 
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I. INTRODUCTION 

It has long been a subject of discussion how to capture 
and describe complex dynamical processes of strongly 
correlated many-body systems that are ubiquitous in na- 
ture. Probably one of the most important problems of 
such kind at present is the one faced in neuroscience. In 
what way biological information processing is encoded 
into spikes of neurons in the brain is the greatest mys- 
tery of present neuroscience. In order to capture subtle 
correlation between activities of neurons, information- 
theoretical techniques have been developed and applied 
to experimental signals [Rieke et al., 1998]. 

On the other hand, we may mention Richardson's cas- 
cade picture of fully-developed fluid turbulence as a suc- 
cessful example of capturing the essence of a strongly 
correlated system, which stated that large eddies grad- 
ually fragment into smaller eddies, and eventually dissi- 
pate due to viscosity. Apart from the difficult problem 
of intermittency, this picture was mathematically formu- 
lated into the Kolmogorov 1941 theory, which gave the 
famous —5/3 energy spectrum [Frisch, 1995]. This suc- 
cess was due to the fact that in fluid turbulence, for- 
tunately, a simple hierarchy in spatial scale (or equiva- 
lently, in temporal / energy scales) is formed as a result 
of complex nonlinear dynamics of a number of modes. 
Since the success in fluid turbulence, attempts have been 
made to capture the dynamics of other kinds of many- 



body systems such as spatio-temporal chaos [Ikeda & 
Matsumoto, 1989], chaotic systems with large degrees of 
freedom [Kaneko, 1994], and stock markets [Arneodo et 
al., 1998a], based on the picture that energy, or more 
generally information or causality, flows from scales to 
scales. 

In analyzing space-scale properties of various kinds of 
signals, the wavelet transform is expected to be a useful 
tool due to its localized nature both in space and scale, 
and many applications have been made [Mallat, 1999]. 
Of course, there are many studies that tried to capture 
the cascade process of fully- developed turbulence directly 
using wavelets [Argoul et al., 1989; Yamada & Ohkitani, 
1991; Arneodo et al., 1998b]. But most of the studies 
treated one-dimensional time sequences of the velocity 
field measured at a single point in fluid, and rather a 
limited number of studies actually investigated the tem- 
poral evolution of correlation between scales of spatially 
extended fields [Toh, 1995]. 

In this paper, we propose a simple visualization 
method of such spatio-temporal correlation between 
scales of spatially extended fields using wavelets, and ap- 
ply it to two typical spatio-temporally chaotic systems. 
As such chaotic systems, we investigate ordinary diffu- 
sively coupled complex Ginzburg-Landau oscillators [Ku- 
ramoto, 1984; Bohr et al., 1998] and somewhat unusual 
non-locally coupled complex Ginzburg-Landau oscilla- 
tors [Kuramoto, 1995]. While the former system exhibits 
spatially smooth amplitude patterns, the latter system 
exhibits fractal spatial patterns. Therefore, the dynami- 
cal processes behind those systems are expected to differ 
from each other significantly, and it is interesting to see 
whether their difference can be detected by our method 
or not. 



II. COUPLED COMPLEX GINZBURG-LANDAU 
OSCILLATORS 

The ordinary diffusively (i.e., locally) coupled complex 
Ginzburg-Landau equation [Kuramoto, 1984; Bohr et at, 
1998] is given by 

W(x,t) = W - (1 + ic 2 )\W\ 2 W + D{\ +ici)V 2 W, (1) 

which can be derived, for example, from equations of os- 
cillatory media in the vicinity of their Hopf bifurcation 
points by the center-manifold reduction technique. Here, 



1 



W(x,t) is a complex amplitude of an oscillator at po- 
sition x and at time t, ci and ci are real parameters, 
and D is a diffusion constant, ft is well known that this 
equation exhibits spatio-temporal chaos in some appro- 
priate parameter region. We call this equation "LCGL 
equation" hereafter. 

In this paper, we treat only spatially one-dimensional 
cases. We fix the length of the system at L = 1, and 
assume a periodic boundary condition. The diffusion 
constant is set at D — 0.0035, and the parameters are 
fixed at C\ = —2 and c 2 = 2. With these values, the 
spatially uniform solution of the LCGL equation is un- 
stable, and the system exhibits spatio-temporal chaos. 
Typical time scale of the system is roughly estimated as 
the time needed for an uncoupled free oscillator to go 
around its limit cycle, and is given by 2ir/c2 — 3. The 
numerical simulation was done in wavenumber space us- 
ing N = 2 10 — 2 14 modes by the pseudo-spectral method 
with a time step of 0.01 (Euler integration). 

Due to the existence of a diffusion term, the solution 
of the LCGL equation necessarily possesses a charac- 
teristic minimal length scale, below which fluctuations 
are strongly depressed. Thus, the amplitude pattern 
t) \ of the solution has a smoothly modulated shape 
as displayed in Fig. |l[ The shortest wavelength deter- 
mined by the diffusion constant is about y/~D ~ 0.06 ~ 
1/16. For the sake of comparison, this value is equated 
with the coupling length of the non-locally coupled sys- 
tem that we explain later. Figure ^| displays temporal 
evolution of the amplitude pattern, and Fig. || shows its 
power spectrum. Since short-wavelength components are 
strongly depressed, the power spectrum decays exponen- 
tially. 
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FIG. 1. Snapshot of the amplitude field W(:r,t) of the 
LCGL equation. N — 2 10 modes are used in the numerical 
calculation. 




FIG. 2. Evolution of the amplitude field of the LCGL equa- 
tion. The horizontal axis indicates the spatial coordinate 
x (0 < x < 1), and the vertical axis the time t (0 < t < 125). 
Time increases upwards. Initial transient is discarded. 
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FIG. 3. Power spectrum of the amplitude field of the LCGL 
equation. The inset shows semi-log plot of the same spectrum. 



2 



As another spatio-temporally chaotic system, we in- 
vestigate a different type of coupled complex Ginzburg- 
Landau oscillators, which was first introduced by Ku- 
ramoto [1995]. Instead of a diffusive interaction, it has 
a non-local interaction; each oscillator feels a non-local 
mean field of other oscillators through a kernel 
which is a decreasing function of Hereafter we use 
<7(|a?|) = go exp(— |a;|/7) as the kernel, where 7 gives the 
coupling range, go is some appropriate normalization 
constant, and assumed to be 27. The equation for this 
system of non-locally coupled oscillators is given by 

W(x,t) = W- (1 + ic 2 )\W\ 2 W 

+K(l + ia) J dx'g(\x' -x\)[W(x',t)-W(x,t)}, (2) 

which we call "NCGL equation" hereafter. This equation 
can be derived, for example, from a model of oscillatory 
biological cells that are interacting through some diffusive 
chemical substance with a finite decay rate. 

Although the difference between this NCGL equation 
and the LCGL equation is only the last interaction term, 
it was shown that this NCGL equation exhibits spatio- 
temporal chaos with remarkably different features from 
the LCGL equation, such as fractal amplitude patterns, 
power-law spatial correlation, and power-law spectrum. 
We again fix the length of the system at L = 1, and as- 
sume a periodic boundary condition. The coupling range 
is fixed at 7 = 1/16, which is equal to the shortest charac- 
teristic length y/~D ~ 1/16 of the LCGL equation. (More 
precisely, by assuming the smoothness of the amplitude 
field W(x' , t) in the interaction term and expanding it up 
to the second order in x' — x around x, we obtain an effec- 
tive diffusion constant D e f f = Kj 2 for the NCGL equa- 
tion. However, the smoothness of the amplitude field is 
not always guaranteed, hence this correspondence is only 
formal.) The parameter values are again set at c\ = —2 
and C2 = 2. With these values, the uniform solution of 
the system is always unstable, and the system behaves 
in a chaotic manner. Further, the behavior of the system 
strongly depends on the remaining coupling strength K; 
the amplitude pattern of the system is smooth for large K 
(~ 1.5), while almost discontinuous for small K (~ 0.5). 
Between these values, there exists a parameter region 
where the amplitude pattern is fractal, and its spatial 
correlation exhibits power-law behavior. We fix K = 0.9 
hereafter, where the system is exactly in this "anoma- 
lous" spatio-temporally chaotic regime. The numerical 
simulation was done in real space using N = 2 10 — 2 14 
oscillators by the 4th-order Runge-Kutta method with a 
time step of 0.05. 

Figure || displays a typical snapshot of the solution 
of the NCGL equation, and Fig. || shows its temporal 
evolution. In contrast to the case of the LCGL equa- 
tion, the amplitude pattern is not completely smooth, 
but is composed of patches of smooth coherent regions 
and strongly disordered regions. The power spectrum 
displayed in Fig. ^ also indicates the difference clearly. It 



exhibits power-law decay rather than the exponential de- 
cay in the case of the LCGL equation, which implies the 
strong anomaly of the amplitude field. Actually, it was 
shown that the amplitude pattern of the NCGL equation 
is fractal, and its fractal dimension varies with the cou- 
pling strength K. These properties are essentially the 
results of the absence of length scale shorter than the 
coupling range 7, and can be explained using a simple 
multiplicative stochastic model to a certain extent [Ku- 
ramoto & Nakao, 1996]. 
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FIG. 4. Snapshot of the amplitude field |W(»,i)| of the 
NCGL equation. N = 2 10 oscillators are used in the numerical 
calculation. 




FIG. 5. Evolution of the amplitude field of the NCGL equa- 
tion. The horizontal axis indicates the spatial coordinate 
x (0 < x < 1), and the vertical axis the time t (0 < t < 125). 
Time increases upwards. Initial transient is discarded. 
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FIG. 6. Power spectrum of the amplitude field of the 
NCGL equation. The inset shows log-log plot of the same 
spectrum. 



Thus, the difference in the interaction term greatly 
changes the amplitude pattern of the system, even 
though the oscillators and characteristic length scales 
are the same. Especially, the power-law behavior of the 
power spectrum reminds us of the —5/3 energy spectrum 
of fluid turbulence. (However, the exponent changes with 
the coupling strength K for our NCGL equation.) The 
powerdaw behavior of the energy spectrum in fluid tur- 
bulence is related to the cascade process of breakdown of 
vortices. Thus, we may naively expect that our NCGL 
equation also possesses some kind of cascade process. 
Of course, the NCGL equation is strongly dissipative 
and there is no such conservative quantity like energy 
or enstrophy. However, we still expect a cascade pro- 
cess of some quantity from long-wavelength components 
to shorter- wavelength components, which may be called 
causality or information, or merely correlation. In the 
following part of this paper, we develop a simple method 
to visualize the spatio-temporal correlation between fluc- 
tuations at different scales, in order to prove the existence 
of such cascade process. 



III. WAVELET-BASED CORRELATION 
ANALYSIS 

Let us consider decomposing the spatial pattern into 
various scales using wavelets, and analyzing the tempo- 
ral correlation between them. We decompose a spatial 
pattern V(x) using an orthonormal wavelet basis as 



V(x) 



a j,fc^,fc(V), 



(3) 



j— — oo k— — oo 



where aj k is an expansion coefficient, and ipj,k(x) is a 
child wavelet generated from the mother wavelet tp(x) by 
translation and dilation as 



where j is a scale parameter, and k is a translation pa- 
rameter. The child wavelets are mutually orthogonal as 



(ipj,k,ipj>,k> 



ipj,k(x')ipj, k ,(x')dx' = 6j,j>6 k<k >, (5) 



and form a complete orthonormal set. Since the child 
wavelet ipj t k(x) is localized both in space and scale 
roughly at position fc/2- 7 and scale 1/2 J , the expansion 
coefficient given by 



a>i,k 



V(x'W k (x')dx' 



(6) 



quantifies the magnitude of fluctuation of V(x) around 
this space-scale point. Hereafter we use Meyer's wavelet. 
It is a real analytic function that decays faster than any 
power function as \x\ — > oo, and its moment of any order 
vanishes, i.e., 



x n ip(x)dx = 0, (n > 0). 



(7) 



Further, it has a smooth Fourier transform with a com- 
pact support. These features are very preferable from 
the viewpoint of physicists. For the details of Meyer's 
wavelet and efficient numerical algorism, see Refs. [Ya- 
mada & Ohkitani, 1991; Mallat, 1999]. Meyer's mother 
wavelet and one of its child wavelet are shown in Fig. ^. 




i> jtk (x) = V&<ip(2 j x-k), 



(4) 



FIG. 7. Meyer's wavelets; mother wavelet (j = 0, k = 0) 
and one of its child wavelet (j = 1, k = 4). 



In the practical numerical calculation, the amplitude 
pattern is a periodic function on [0, 1) represented by 
N = 2™ discrete points. Therefore, we extend the pattern 
over the whole real axis by a periodic extrapolation in 
order to apply the wavelet transform. Consequently, in 
the actual numerical calculation, the first subscript j in 
the summation of Eq. (|^) indicating the scale runs from 
to n — 1, and the second subscript fc indicating the 
translation runs from to 2 J ' — 1. 

Rather than using the obtained expansion coefficient 
directly, we define a new coarse-grained field bj(x) (0 < 
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x < 1) by a piecewise constant interpolation from the 
expansion coefficient a,j t k as 



bj(x) = log|oj,fc| 



<x< — — , fc = 0,l,...,2 3 - 1 



23 



(8) 

Here we take a square of the coefficient, because we are 
interested in the absolute intensity of the fluctuation at a 
certain position and scale. Further taking the logarithm 
is merely for convenience sake here; we obtain similar re- 
sults without taking the logarithm. However, there also 
exists some discussion claiming that taking the logarithm 
is more appropriate when the system under consideration 
possesses a random multiplicative process [Arneodo et 
at, f998a; Arneodo et al, 1998b], and the NCGL equa- 
tion is actually considered this case. Figures || and || 
display the coarse-grained fields bj(x) for the LCGL and 
NCGL equations. The obtained data for the LCGL equa- 
tion are very noisy for large j values, since components 
whose wavelengths are shorter than the dissipation scale 
VB ~ 1/16 decay very quickly. On the other hand, we 
can observe clear positional correlation between scales 
down to a very small wavelength in the case of the NCGL 
equation. 




FIG. 8. Coarse-grained field bj(x) obtained from the loga- 
rithm of the squared wavelet coefficient of the LCGL equation. 
The top panel shows the original amplitude pattern |W(a?)|. 
The scale j ranges from 1 to 9 from above. Each graph is 
normalized with its mean and variance. N = 2 10 modes are 
used in the numerical calculation. 
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FIG. 9. Coarse-grained field bj(x) obtained from the loga- 
rithm of the squared wavelet coefficient of the NCGL equa- 
tion. The top panel shows the original amplitude pattern 
|W(a;)|. The scale j ranges from 1 to 9 from above. Each 
graph is normalized with its mean and variance. N = 2 10 
oscillators are used in the numerical calculation. 



Let V\{x) be the amplitude pattern at time t, i.e., 
V\(x) — \W{x,t)\, and V^(x) be the amplitude pattern 
evolved from Vx(x) for a period of At, i.e., V%(x) = 
|W(a;,i+At)|. From these amplitude patterns, we obtain 
coarse-grained fields 6^ (x) and 6| 2 (x) at scales ji and j2- 
We then define deviations of these fields from their mean 
values as 



A6 J \(x) = 6j l (x)-(6) i ,7), 



and 



Abl( x )=bl(x) 



where I is a constant-valued function I(x) = 
inner product of /(x) and g(x) is defined as 



(f,g) = / f(x)g(x)dx. 
Jo 



(9) 

(10) 
1, and the 

(11) 



Now we define 



kind of cross-correlation matrix 



Cj 1 j 2 (At) from these A6j i (x) and Airj 2 (x) as 



(Ab^, A6jJ (Abl, Abl) 



1/2- 



(12) 
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From this definition, Cj 1 .j 2 (At) — C j2jl (— At) obvi- 
ously follows. The overlines in the above equations indi- 
cate temporal average, assuming the stationarity of the 
spatio-temporal chaos under consideration. 

When At = 0, V?(x) represents the same amplitude 
pattern as V\{x). Thus Abj 1 (x) and Ab 2 j 2 (x) are iden- 
tical, and Cj 1 .j 2 is merely a symmetric matrix. When 
At 7^ 0, V2(x) differs from Vi(x) to some degree depend- 
ing on At. Therefore Cj lj2 is no longer symmetric, and 
its asymmetry is expected to characterize the variation of 
spatio-temporal correlation between scales. In order to 
see this asymmetry clearly, let us decompose the matrix 
Cj lt j 2 into symmetric and antisymmetric parts as 

D jl>j2 m = \ {C juh (At) + C htjl (At)} , (13) 

and 

E juh (At) = C juh (At) - C h , h (At). (14) 

The symmetric part Dj 1 j 2 is invariant under time- 
reversal transform At <-» —At, while the antisymmet- 
ric part Ej lt j 2 is not invariant (antisymmetric). Thus 
the matrix Ej 1 j 2 is expected to quantify the variation 
of spatio-temporal correlation between scales that is not 
symmetric to time-reversal, namely, to give a certain 
measure of the flow of information or correlation between 
scales. 

In Fig. |l0|, the symmetric part Dj lt j 2 and the antisym- 
metric part Ej 1 _j 2 of the correlation matrix obtained for 
the LCGL equation are displayed using color codes for 
several values of the time difference At. The top row 
displays results obtained for the same time (At = 0), 
and the lower rows display results for larger values of the 
time difference At. The vertical axis of each graph rep- 
resents the scale j\ of the pattern at the reference time, 
and the horizontal axis represents the scale j'2 of the pat- 
tern At after the reference time. First, for At = 0, the 
antisymmetric part is merely zero and only the symmet- 
ric part takes positive values. The diagonal components 
naturally take largest values, since the self-correlation of 
the fluctuation at the same scale is strongest. Besides, 
it can be seen that there exists rather strong correlation 
between adjacent scales in the dissipation region deter- 
mined by the diffusion constant D (the region with the 
scale slightly smaller than j ~ 4). When the time dif- 
ference At becomes a little larger, the symmetric part 
diminishes. But still the diagonal components and the 
dissipation region maintain higher correlations than the 
others. Now, if we look at the antisymmetric part, there 
appear small localized regions in the neighborhood of the 
dissipation region across the diagonal, which have nega- 
tive and positive correlations, respectively. This implies 
that fluctuation at a certain scale at the reference time 
(vertical axis) possesses relatively strong correlation to 
fluctuation at slightly smaller scale At after the reference 
time (horizontal axis). In this case, it can be considered 
as visualizing the dissipation due to the diffusion term. 



As the time difference At becomes further large, both the 
symmetric and antisymmetric parts tend to take smaller 
values, and correlation between the patterns vanishes. 
The intensity of the antisymmetric part is typically 20% 
of the symmetric part at its maximum. 

Figure [ll] displays the symmetric part Dj 1 J2 and the 
antisymmetric part Ej lt j 2 of the correlation matrix ob- 
tained for the NCGL equation for several values of the 
time difference At using color codes, as in the case of 
Fig. |l(| When there is no time difference (At = 0), 
only the symmetric part takes finite values, and its di- 
agonal components take largest values as in the case of 
the LCGL equation. However, the region where the cor- 
relation between adjacent scales takes its peak value is 
located at much smaller scale. When the time differ- 
ence At becomes a little larger, the difference becomes 
more remarkable; there appear wide coherent regions 
with positive and negative correlation, which spread over 
scales shorter than the characteristic length of the system 
(J ~ 4) determined by the coupling range 7. This asym- 
metry is maintained up to a considerably large value of 
At, which indicates that long-wavelength components at 
the reference time maintain strong correlations widely 
to shorter-wavelength components at later times. As 
we increase the time difference At further, the corre- 
lation between scales gradually vanishes from the long- 
wavelength region. But it takes much longer than the 
case of the LCGL equation. These facts suggest the exis- 
tence of cascade-like propagation of some quantity from 
long wavelength to shorter wavelength. The maximum 
intensity of the antisymmetric part is about 10% of the 
symmetric part. 



G 





FIG. 10. Evolution ol the symmetric part Dj 1: j 2 (left col- 
umn) and the antisymmetric part Ej 1 j 2 (right column) of the 
correlation matrix with the time difference At for the LCGL 
equation. In each figure, the vertical axis indicates the scale 
ji at the reference time, and the horizontal axis indicates the 
scale j2 at the time At after the reference time. The scales 
of both axes range from to 13. Time intervals between pat- 
terns are At = 0, 1, 2, 4, 8 from the top. TV = 2 14 modes are 
used in the numerical calculation. The color scale of Dj lt j 2 
is from (blue) to 1 (red), while that of Ej ± j 2 is from —0.1 
(blue) to 0.1 (red). 
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FIG. 11. Evolution of the symmetric part Dj lt j 2 (left col- 
umn) and the antisymmetric part Ej 1 j 2 (right column) of the 
correlation matrix with the time difference At for the NCGL 
equation. In each figure, the vertical axis indicates the scale 
ji at the reference time, and the horizontal axis indicates the 
scale j2 at the time At after the reference time. The scales of 
both axes range from to 13. Time intervals between patterns 
are At = 0, 2, 5, 20, 40 from the top. N = 2 14 oscillators are 
used in the numerical calculation. The color scale of Dj lt j 2 
is from (blue) to 1 (red), while that of Ej lt j 2 is from —0.05 
(blue) to 0.05 (red). 
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In order to see the temporal evolution of the antisym- 
metric part of the correlation matrix in more detail, an- 
tisymmetric components Ejj+i between adjacent scales 
j and j + 1 normalized by the maximum value at each 
scale are displayed in Figs, [l^ and [l3| for scales smaller 
than the characteristic length. In the case of the LCGL 
equation, there is no clear ordering of the temporal vari- 
ation of correlation. Each curve has its maximum at 
some small value of At, and vanishes quickly. In the case 
of the NCGL equation, on the other hand, there exists 
a clear time ordering of the curves in scale; the correla- 
tion between long- wavelength fluctuation (small j values) 
becomes large earlier and decreases quickly, while the 
correlation between short-wavelength fluctuation (large 
j values) becomes large later and decreases very slowly. 
(Note the difference of time scale between the graphs for 
the LCGL and NCGL equations.) 
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FIG. 12. Antisymmetric part of the correlation coefficient 
Ejj+i between two successive scales j and j + 1 as a function 
of the time interval At for the LCGL equation. Each curve is 
normalized by its maximum value. 
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FIG. 13. Antisymmetric part of the correlation coefficient 
Ejj+i between two successive scales j and j + 1 asa function 
of the time interval At for the NCGL equation. Each curve 
is normalized by its maximum value. 



IV. DISCUSSION 

As we demonstrated, our simple method based on the 
correlation matrix seems to visualize the spatio-temporal 
correlation between scales of spatio-temporal chaos suc- 
cessfully. Our method clearly visualized the difference 
of dynamical process between two systems with different 
natures. Especially, for the non-locally coupled system, 
it verified the existence of cascade-like propagation of 
correlation from long-wavelength components to shorter- 
wavelength components. 

First, the use of the wavelet transform is essential for 
our results. It seems difficult to obtain similar results 
using a method that flattens out the spatial information 
completely, like the Fourier transform. As can clearly 
be seen from the definition of the inner product given 
in Eq. dTl|), our method captures the simultaneous fluc- 
tuation of two fields at the same position in space. For 
example, if we use a different definition of the inner prod- 
uct that loses the spatial information by taking spatial 
average of each field first, we cannot observe the cor- 
relation between scales so clearly as in Figs. [lO] and [n]. 
Similarly, if we change the definition of the inner product 
so as to multiply the fields with their origins shifted by 
half the system size, the correlation matrix almost van- 
ishes. (This tendency is more remarkable for the NCGL 
equation. In the case of the LCGL equation, though not 
so clear, we obtain similar asymmetric correlation matri- 
ces with these modified definitions of the inner product. 
This is because with the values of D and L we used in our 
calculation, the ratio of the diffusion constant to the sys- 
tem size is considerably large, and some coherence over 
the whole system still remains. If we further decrease the 
value of -D, the correlation matrix tends to vanish as in 
the case of the NCGL equation.) 

Conversely, the fact that we succeeded in visualization 
using the definition of the inner product given in Eq. ( |l"l| ) 
implies that the cascade process occurs in spatially local- 
ized regions in our system. However, it is frequently seen 
in other systems that the spatial structure of the pattern 
collapses while moving constantly, as in the case of con- 
vective instability. In order to detect such phenomena, 
it would be possible to generalize the definition of the 
inner product in such a way that the distance between 
the origins of the fields is increased with At. 

Precisely speaking, the result visualized by our method 
is not a flow of causality, but merely a flow of correlation. 
Namely, it indicates that fluctuation at the reference time 
at some scale varies in unison with fluctuation at the 
later time at some other scale, but it does not necessar- 
ily mean that the former one actually affects the latter 
one. Of course, it is natural to interpret it as causality 
in our cases. But in order to be exact, it will be possible 
to perturb a certain mode at some scale, e.g. by using a 
sinusoidal wave, and visualize the propagation of its af- 
tereffect to prove that it is indeed a causal relationship. 

Finally, though we did not give a detailed discussion in 
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this paper, we will be able to investigate the dynamical 
process of the spatio-temporal chaos of non-locally cou- 
pled oscillators in more detail by quantitatively analyzing 
the results obtained by our method. For example, from 
the peak positions of the curves shown in Fig.|T5|, 

we will be able to study the dependence of characteristic 
time scale of the underlying cascade process on the cou- 
pling strength. Such a study will also be an interesting 
future subject. 
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